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Recently interesting data on dNcYi/dr] in Au-Au collisions [r] = — lntan(0/2)) with 
the centrality cuts have been reported by PHOBOS Collaboration. Their data are 
usually divided by the number of participants (nucleons) in collisions. Instead of 
this way, using the total multiplicity A'ch = /(diVch/d7?)d?7, we find that there is 
scaling phenomenon among {N dN ch_ / drj = dn/drj with different centrality cuts 
at -\/snn = 130 GeV. To explain this scaling behavior of dn/dr], we consider the 
stochastic approach called the Ornstein-Uhlenbeck process with two sources. More- 
over, comparisons of dn/dr/ at ^snn = 130 GeV with that at -^/snn = 200 GeV 
have been made. A possible detection method of the quark-gluon plasma (QGP) 
thorough dn/dr? is presented. 



1 Introduction 



Recently PHOBOS Collaboration has publislied an interesting data on dNch/dr] 



ij] = — lntan(6'/2))Q in Au-Au collisions at ^snn = 130 GeV [1]. The authors 
of ref. [1] have calculated the following quantity, 
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^^"'^ ■ /(A^part, Ncolh V) 



{Np^rt)/2 dr] 

where A^part and A'cou mean the number of participants (nucleons) and number 
of colhsion particles in Au-Au colhsions. It depends on the centrahty cuts. The 
function /(A^part, A^'cou, r] = 0) is an increasing function, as A'part increases. 

In this paper, instead of eq. (1), we consider the following physical quantity. 



1 dA^ch _ dn 

A^ch dr] ~ df]' ^ ^ 

where A'ch = J{dNch/dri)dri, and J{dn/dri)dri = 1. In Fig. 1, two sets of dn/dr] 
are shown. They suggest us that there is scaling among dn/dr^'s with different 
centrality cuts. Thus dn/dr] is named a kind of the probability density, because 
dn/dr] = f{r]), where dn/dr] is function of r] only This fact probably implies 
that the stochastic approach is available in analyses of dn/dr? pP] 

Contents of the present paper are as follows. In the second paragraph, the 
Fokker-Planck equation is considered. In the third one, concrete analyses are 
presented. In the final one, concluding remarks are given. 



^ Multiplicity distributions in high energy collisions, i.e., P{n, (n))'s are the prob- 
ability distributions, which are function of n and (n). It is known that the KNO 
scaling functions, 

lim {n)P{n) = 'il^{z = n/ {n)) 

are described by solutions of various Fokker-Planck equations [2-4]. 

^ Dokshitzer has calculated the generalized gamma distribution in QCD [5] 



P{z = n/{n))^^^-^^^e^p[-{Dzr] , 

where z, /i, D, 7 are KNO scaling variable, 1 — 7 = 1//U, a parameter, the anomalous 
dimension in QCD, respectively. This is a steady solution of the following Fokker- 
Planck Equation 
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Fig. 1. Two sets of dn/dr] with different centrality cuts [1]. Data with 15-25% are 
used in both Figures, for the sake of comparisons. 

2 Stochastic approach to dn/dr] 



It is well known that the rapidity {y ps ry) is a kind of the velocity. Moreover, 
in collisions leading particles in the beam and target nuclei, i.e., nucleons in 
gold at RHIC experiments, collide each other, lose their energies and emit 
various particles. In other words, there are fluctuations in the rapidity space 
or the pseudo-rapidity {r}). Thus we would like to assume that a damping law 
governs the rapidity space 



^ = -7^ + Mt) , (3) 

where X — y or rj, ^ and /w(0 are the frictional coefficient and the white 
noise, respectively. Prom eq. (3), we can derive the following Fokker-Planck 
equation. 



dP{x, t) 
dt 



= 7 



d_ 1^2^ 

dx 27 dx'^ 



P{x, t) , 



(4) 



where <T and P{x, t) are the dispersion and the probability density, respec- 
tively. Equation (4) is called the Ornstein-Uhlcnbeck (0-U) process [6,7]. To 
look for the solutions of eq. (4), we have to assume the source functions at 
"time" t — Q. The typical solution with 5(77 — r/o) = P{r], t = 0) is given as 



3 



exp 



(t) - rjoe 



(5) 



where V^lt) = {a'^/2j){l — e ^'^*) and r/o is the initial rapidity. 



In Fig. 2(a), wc depicted a simphfied picture of heavy-ion colUsions. According 
to Fig. 2(a), a model of two sources 0.5 x 5{rj — r/max) and 0.5 x 6{ri + ?7max) at 
t = seems to be reasonable. As our present model is very simple, 0.5 x Nch 
particles have been produced at r/max and the same 0.5 x A^ch particles at — r/max 
at t — 0. In this case we have the following expression which is the sum of two 
solutions. 
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(6) 



The evolution of eq. (6) is shown in Fig. 2(b). 
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Fig. 2. (a) Simplified picture for A-A collisions. Tlie black circles mean nucleons. 
(b) Evolution of eq. (6) with two sources at r^max and — r/max- 



3 Analyses of dn/dr^ by means of eq. (6) 



By making use of eq. (6), we can analyse dn/drj shown in Fig. 1. Our results 
are shown in Fig. 3 and Table 1. In Fig. 4, we examine whether or not the 
dispersion ¥"^{1) and p — 1 — e~'^'^ depend on the centrahty cuts. As is seen 
in Figs. 3 and 4, the scahng behavior among dn/dr^'s at ^snn = 130 GeV is 
explained by eq. (6) with small changes in the dispersion V'^{t). It can be said 
that the scaling behavior is explained by the 0-U process with two sources 
at the beam (j/b or r/max) and target (j/t or — ?7max) rapidities. Of course, it is 
obvious that the single source cannot explain it. 
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Fig. 3. Analyses of dn/dr] by eq. (6). See Table 1. 
4 Concluding remarks 

First of all, it can be stressed that there is the scaling among d?2/d?7's with 
various centrality cuts at ^snn = 130 GeV. 

Second the scaling behavior of dn/drj is described by the solutions of the 
Fokker-Planck equation, i.e., eq. (6)[^. This suggests that dn/drj with the cen- 
trality cut 0-6% do not show singular /or particular phenomenon relating to 
signatures of the Quark-Gluon Plasma (QGP). Of course, we should pay our 
attention that we are handling the averaged quantity in statistics. At present, 
however, it is difficult to conclude that the QGP is created, and the signature 

^ To estimate the "thermalization time" of the QGP, Hwa has considered the 
Fokker-Planck equation for the motion of the quarks and gluons in nuclei [8]. See 
also ref. [9] in which the Wiener process is considered for the quarks and gluons. 
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Table 1 

Estimated parameters in our analyses by eq. (6) with two somxes. Evolution of eq. 
(6) is stopped at minimum x^'s- = 0.006 ~ 0.004. 



Fig. 3 
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centrality (%) 


45-55 


35-45 


25-35 


15-25 


6-15 


0-6 


P 


0.872±(5p 


0.875±5p 


0.878±5p 


0.882±(^p 


0.886±5p 


0.888±(5p 


vHt) 


3.83±0.27 


3.61±0.21 


3.23±0.16 


3.00±0.13 


2.72±0.10 


2.47±0.08 


(A^part) 




93 


135 


197 


270 


340 


iVch 


662±10 


1056±16 


1582±23 


2270±34 


3199±49 


4070±63 


X^/n.d.f. 


8.61/51 


7.63/51 


5.88/51 


5.35/51 


3.57/51 


3.82/51 
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Fig. 4. Dispersion V'^{t) and p of Fig. 3 and Table 1. 

from the QGP are washed out by the strong interactions between hadrons, if 
the QGP is created. 

Here we should carefully observe Fig. 1. As is seen in Fig. 1, there are small 
differences in dn/dr^'s over the range |?7|^2 with centrality cut 0-6% and oth- 
ers. To explore the differences more carefully, we need dn/dr] with smaller 
centrality cuts, 0-3% ~ 0-5%. 

Moreover, we can add the following fact. Very recently, PHOBOS Collabora- 
tion has reported the data on dNch_/drj with centrality cut 0-6% at ^snn = 200 
GeV [10]. They are compared with dn/dr] at i/snn = 130 GeV in Fig. 5. 
Roughly speaking, the scaling on dn/drj holds between ^snn = 130 GeV and 
200 GeVQ. The distribution of dA'ch/d?] or dn/d?7 with the centrality cut 0- 
6% does not show the particular behavior relating to the QGP in the sense of 

^ Using the Bjorken's picture [11] for the calculation of energy density near |Ar/| < 
0.5 with the geometrical picture of the gold {R-r « 6-7 fm, ctq « 1-2 fm, V ~ 
7ri?|(cro) ^ 300 fm^), we obtain the following values 
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Fig. 5. Comparisons of dn/dr) with the centrahty cut 0-6% at y^SNN = 130 
GeV and 200 GcV. Solid hne is obtained for latter energy: p = 0.879 ± 0.007, 
= 2.67 ± 0.24, x^/n-d.f = 0.63/22. 

average. To investigate the particular phenomena like the turbulence and/or 
deflagration in dA^^ch/dr/, we need to analyse the single event with smaller 
centrality cut than 0-6%. 

Moreover, analyses of event-by-event, for example the intermittency and the 
wavelet, seem to be necessary for dNch/dr] with smaller centrality cut [12-18]. 
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